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Figure 1: Types of Feynman graphs encountered in a tree-level calculation of weak boson pair production in 
association with two jets. Wavy lines depict weak bosons and spiralled lines are gluons. Quark flavour labels 
and weak boson labels (W + /W~ /Z°) have been deliberately suppressed to highlight the general topology of 
the graphs which may contribute to different processes. 



1. Introduction and Motivation 

Precision calculations of standard model processes are essential for interpreting the signals 
measured at the Large Hadron Collider (LHC) and for fully realising the discovery potential of 
this experiment. Next-to-leading order (NLO) calculations in perturbative QCD have proved very 
successful when used in Tevatron analyses, and are a good way to reduce theoretical uncertainties 
in a description of a given process. On top of this, merging an NLO calculation with a parton 
shower provides a realistic hadron-level prediction for an event whilst maintaining NLO accuracy 
for inclusive observables. 

In this talk I shall discuss the production of a pair of weak bosons in association with jets - 
specifically the two processes pp — > W + W + jj and pp — > W + W~ jj. I will describe the compu- 
tation of the NLO QCD corrections to both processes, as well as the merging of pp — > W + W + jj 
with a parton shower, done in the framework of the POWHEG BOX [Q]. 

Figure |l| depicts the structure of some of the types of Feynman graphs one encounters in a 
tree-level calculation of weak boson pair production along with two jets. All of these types of 
graph contribute towards the process pp — > W + W~ jj, but only graphs of the type shown in the far 
right of the figure contribute towards pp — > W + W + jj. Here, charge conservation requires the two 
W + bosons to be emitted from separate quark lines and this leads to an unusual theoretical property 
- the cross section for this process remains finite even if the requirement that two jets are observed 
is lifted. This will be investigated later on, and I will present results for W + VK + + n jets, where 
n = 0, 1,2. The calculation of pp -4 W + W + jj can be seen as a stepping stone to the calculation of 
pp — > W + W~jj, since it involves a small subset of the Feynman graphs needed for the latter. 

Both processes are 2 — ^ 4 processes, and to calculate the QCD corrections to them one needs 
to deal with one-loop, six-point tensor integrals of relatively high rank. There is thus a theoretical 
incentive in performing these calculations and much progress has been made over the past few 
years in the methods used to compute them - this will be discussed in the following section. But 
before this, I will go on to discuss the study of both processes at the LHC in a bit more detail. 

1.1 W + W + jj at the LHC 

At y/s =14 TeV, the cross-section for this process is about 1 pb (40% of this for W~W~jj) and 
is therefore accessible. In the following we take the W + bosons to both decay leptonically, giving 
rise to a nearly background-free signature which involves same-sign leptons. This is an interesting 
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process to study in its own right, but there are other reasons to study it: pp — > W + W + jj is a back- 
ground to physics both within and beyond the standard model. For example, it is possible to use 
same-sign lepton pairs to study double parton scattering at the LHC to which pp — > W + W + jj 
is a background. Beyond the standard model, resonant slepton production in R-parity violating 
SUSY models [g], diquark production [fjj, and doubly charged Higgs boson production [Jsj] are ex- 
amples of processes which also lead to a signature of two same-sign leptons, missing energy, and 
jets. 

1.2 W + W~jj at the LHC 

The production of a boson pair in association with zero, one or two jets is an im- 

portant background to Higgs boson production, especially when the decay H — > W + W~ opens up. 
Although most of the sensitivity in Higgs boson searches comes from the zero jet processes, which 
have the largest cross-section, the production of a Higgs boson in association with two jets is also 
relevant - about 10% of Higgs events at the LHC involve two jets [^, |XC| ]- The production of a 
Higgs boson via weak boson fusion (WBF) also has a sizeable cross-section. The signature of this 
process includes two forward tagging jets and pp — > W + W~jj is an irreducible background to this. 
As we did for W + W + jj, in the following we will take both W bosons to decay leptonically. The 
resulting signature of two opposite-sign leptons, jets and missing energy is also a background to a 
classic beyond the standard model physics search. 



2. Method of calculation 

2.1 The NLO QCD corrections 

NLO QCD calculations of processes involving more than five particles is difficult. For the 
virtual amplitude, the number of Feynman diagrams needing evaluation grows factorially with the 
number of particles in the process. In addition to this, the one-loop tensor integrals which need to 
be computed become more involved. However, a refinement of traditional computation methods, 
as well as the development of new techniques based on unitarity and on-shell methods, have seen 
a significant growth in the number of 2 — > 4 processes (and even a 2 — > 5 process) known at NLO 



in the past few years (see [11] for a recent review). Platforms for the automation of NLO-accurate 



processes are currently being developed (see e.g. JT2| , |l3| , 14, 15]). 

As described in detail in the papers [jl|, ||], the technique of D-dimensional generalised unitar- 



ity [16] was used to obtain the virtual part of the amplitude for the QCD processes pp — > W + W + jj 
and pp — > W + W~ jj. It is worth pointing out that, as currently formulated, on-shell methods re- 
quire working with an ordering of external lines - these are colour ordered or primitive amplitudes. 
It is only colour-charged particles which are ordered in primitive amplitudes and so all possible in- 
sertions of the colourless weak bosons must be considered for any tree-level or one-loop primitive 
amplitude. The D-dimensional unitarity cuts reduce one-loop primitive amplitudes to products of 
tree-level helicity amplitudes, and a certain amount of difficulty exists in ensuring no over-counting 
takes place when combining the cuts of different parent diagrams. Nevertheless, this is just book- 
keeping and these two calculations demonstrated that unitarity methods can deal with more com- 
plicated, colourless final states. The tree-level helicity amplitudes themselves are calculated using 



Berend-Giele recursion relations [17]. 
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Figure 2: The dependence on factorisation and renormalisation scales of cross -sections for 
pp — > e + fi + V e Vji +n jets, n = 0,1,2 at leading and next-to-leading order in perturbative QCD. Here 



2.2 Merging with a parton shower 

Methods which include both the benefits of an NLO calculation and a parton shower model 
(NLO+PS generators) have become available in recent years - two frameworks are currently being 



used for collider physics: MC0NLO Q18Q and POWHEG Q19p . A general computer framework for 
building a POWHEG implementation of an arbitrary NLO process exists - the POWHEG BOX 
Here, one needs only to supply a few ingredients: phase-space and flavour information, the Born 
and real matrix elements and the virtual matrix elements for a given NLO process. The implemen- 
tation of pp ->■ W + W + jj in the POWHEG BOX is reported in This was the first time a 2 -> 4 
process was implemented in a NLO+PS generator. 

Since all of the ingredients needed by the POWHEG BOX were already known from [JTJ , the 
POWHEG implementation of this process did not present any special problem, except for a non- 
trivial issue of high computational demands coming from the virtual corrections. The technical 
details of how this problem was dealt with are described in detail in [ph. The resulting code is 



public and is available at the website [20]. 



3. Results 

3.1 Selected results for pp -> W + W + jj at the LHC 

First I will present results from the NLO calculation of pp — > W + W + jj, taken from We 
consider proton-proton collisions at a center-of-mass energy y/s = 14 TeV. We require leptonic 
decays of the W -bosons and consider the final state e + pi + v e v ll . The W -bosons are on the mass- 
shell and we neglect quark flavour mixing. We impose standard cuts on lepton transverse momenta 
p±j > 20 GeV, missing transverse momentum /?j_, m i ss > 30 GeV and charged lepton rapidity |t]/| < 
2.4. We define jets using anti-&^ algorithm, with R = 0.4 and with a transverse momentum cut 
p±j = 30 GeV on the two jets. The mass of the W -boson is taken to be mw = 80.419 GeV, 
the width IV = 2.140 GeV. W couplings to fermions are obtained from CKqed^z) = 1/128.802 
and sin 2 dw = 0.2222. We use MSTW08LO parton distribution functions for leading order and 
MSTW08NLO for next-to-leading order computations, corresponding to a s (Mz) = 0.13939 and 
OC s (Mz) = 0.12018 respectively. We do not impose lepton isolation cuts. All results discussed 
below apply to the QCD production pp — > W + W + jj; the electroweak contribution to this process 
is ignored. 
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Figure 3: The kinematic distribution for the transverse momentum of the third hardest jet in the QCD 
production of pp -t e + jx + v e v ll + 2 jets. The pure NLO result and the result with POWHEG+P YTHIA are 
both shown. 



Figure || shows the dependence of the production cross-sections for pp — > e + pl + v e v ll + n jets 
on the renormalisation and factorisation scales, which we set equal to each other. Considering 
the range of scales 50 GeV < fX < 400 GeV, we find the two-jet inclusive cross-section to be 
a LO = 2.7 ± 1.0 fb at leading order and a NLO = 2.44 ± 0. 18 fb at next-to-leading order. The forty 
percent scale uncertainty at leading order is reduced to less than ten percent at NLO. We observe 
similar stabilization of the scale dependence for the 0- and 1-jet exclusive multiplicities. Combining 
these cross-sections we obtain a total NLO cross-section of about 2.90 fb for pp — > e + pL + v e v l i 
inclusive production. This implies about 60 e + pL + + e + e + + events per year at the LHC 

with 10 fb -1 annual luminosity. While this is not a gigantic number, such events will have a very 
distinct signature, so they will definitely be seen and it will be possible to study them. 

The dramatic change in the two-jet exclusive cross-section apparent from figure |2| is discussed 
and investigated in [|l|]. We find that the feature observed here, that the two-jet exclusive is sig- 
nificantly smaller than the two-jet inclusive, remains present when we increase the jet cut and so 
allow for greater perturbative convergence of the exclusive cross section. This smallness implies 
that quite a large fraction of events in pp — > e + n + v e v jl + > 2 jets have a relatively hard third jet. 
This feature may be useful for rejecting contributions of pp — > W + W + jj when looking for multiple 
parton scattering. 

Next I present results from the POWHEG implementation of pp — > W + W + jj, taken from the 
paper Here the set-up is as described above, but we consider pp collisions at a different centre 
of mass energy: y/s = 7 TeV. A dynamic scale is used for the renormalisation and factorisation 
scales: 

Hr = Hf = {p±,i+P±.2+E± ;Wl +E ±}W2 )/2, Ex,w = \l m w + P\,w 

where p±.Wi> P±.w 2 > P±,i an d P±.2 are the transverse momenta of the two Ws and the two emitted 
partons in the underlying Born configuration. 

With no jet cuts, but with the leptonic cuts described above, we find the cross-section for to be 
1.11 ±0.01 fb for the pure NLO result, and a slightly lower cross-section of 1.06 ±0.01 fb when 
events are generated by POWHEG and are subsequently showered with P YTHIA. A comparison 
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Figure4: Left pane: the production cross-section of the process pp — > (W + —> V e e + )(W~ — > fl^V^jj at the 
7 TeV LHC in dependence on the factorisation and renormalisation scales jif = Hr = ji at both LO and NLO 
in perturbative QCD. Right pane: the dependence of the cross-section on centre of mass energy s/s with LO 
results in dashed blue and NLO results in solid red. Three choices of pL are shown: /! = m<w,2m w ,4mw. 



of kinematic distributions was carried out in ^ and for the most part, there was good agreement 
between the NLO and the POWHEG+PYTHIA results. However, there were some distributions 
which showed expected and marked changes, one of which I shall highlight in this talk. Figure || 
shows the transverse momentum of the third-hardest jet. Since at NLO it is only the real radiation 
which contributes to this distribution, we see a divergence for small px,j 3 m the pure NLO result. 
In contrast one can see the Sudakov peak in the POWHEG + PYTHI A result, and the distribution 
goes to zero as pxj 3 0. 

3.2 Selected results for pp -»■ W + W~jj at the LHC 

Here I will present selected results from the calculation of pp — > W + W~ jj, taken from the 
paper [gj. Here the W bosons decay leptonically: W + W~jj — > e + jJ,^v e v^jj. The full results with 
generic opposite-sign leptons can be obtained from these by multiplying by a factor of four. We use 
the same leptonic cuts and electroweak input parameters as were described in the pp — > W + W + jj 
results section above. However, here of course a jet cut must be applied and two jets observed in 
order to obtain a finite cross-section: we take pxj > 30 GeV and |7jy| < 3.2. 

Figure || shows the dependence of the production cross-section on renormalisation and fac- 
torisation scales, which are again set equal to each other, at a centre of mass energy ^/s = 7 TeV. 
The dependence of the cross-section on centre of mass energy is also shown in figure g. One ob- 
serves a dramatic reduction in scale dependence in going from leading order to next-to-leading 
order. Considering a range of scales my? < }X < Amy/ we obtain a cross section at leading order 
Olo = 46 ± 13 fb and at NLO Gnlo = 42 ± 1 fb. Assuming fifty percent efficiency, with 5 fb~ ! of 
data at the 7 TeV run of the LHC, we expect about 400 dilepton events e + ji~ ,e + e~ e~ jx~ . 
It is interesting that at NLO, the dependence of the cross-section on centre of mass energy is 
almost linear. If one defines an 'optimal' scale choice to be the choice of scale for which NLO 
corrections are smallest then this 'optimal' scale shifts from 2mw at 7 TeV to Amy? at 14 TeV. 

Finally I present two kinematic distributions for this process which are relevant for a Higgs 
boson search at the LHC. The left pane of figure |5| plots the relative azimuthal angle between the 
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Figure 5: Distributions of lepton opening angle and jet pseudorapidity difference for the process 
pp — !• (W + -» V e e + )(W~ — > Ll~Vn)jj at the 7 TeV LHC. LO results are shown in blue, NLO results in 
red. The uncertainty bands are for scale mw < fi < 4mw and the solid lines show the results at n — 2mw- 

leptons which peaks at <j) e -n+ = %. This is in contrast to leptons produced via the mechanism 
H — > WW — > e + /i~w where this angle tends to be small. The pseudorapidity difference between 
the two leading jets, AtJji/2 = f]j\ — f]j2, is plotted in the right pane of figure ||. This is a useful 
distribution for studying Higgs boson production via WBF - this mechanism leads to jets which 
tend to have a large |Atj ; i ; 2|- For a Higgs produced via gluon fusion and, as we see here for 
pp — > W + W~ jj, this distribution is peaked around |Atj ; i ; 2| = 0. The significant reduction in the- 
oretical scale uncertainties can also be seen in these distributions, and there is no observed shape 
change in going from LO to NLO. These observations were typical of all kinematic distributions 
considered in [H]. 

4. Conclusion 

In this talk I have presented the NLO QCD corrections for the process pp -^-W + W + jj and 
the process pp — > W + W~jj which were computed using the method of D-dimensional generalised 
unitarity. A significant reduction in the theoretical uncertainties of an LHC prediction is observed 
for both processes. The process pp — > W + W + jj has been implemented in the POWHEG BOX 
which matches the NLO result with a parton shower. I look forward to measurements of pairs of 
weak bosons and jets at the LHC. 
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